Quasi-equilibrium optical nonlinearities in spin-polarized GaAs 
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Semiconductor Bloch equations, which microscopically describe the dynamics of a Coulomb in- 
teracting, spin-unpolarized electron-hole plasma, can be solved in two limits: the coherent and the 
quasi-equilibrium regime. These equations have been recently extended to include the spin degree of 
freedom, and used to explain spin dynamics in the coherent regime. In the quasi-equilibrium limit, 
one solves the Bethe-Salpeter equation in a two-band model to describe how optical absorption is 
affected by Coulomb interactions within a spin-unpolarized plasma of arbitrary density. In this 
work, we modified the solution of the Bethe-Salpeter equation to include spin-polarization and light 
holes in a three-band model, which allowed us to account for spin-polarized versions of many-body 
effects in absorption. The calculated absorption reproduced the spin-dependent, density-dependent 
and spectral trends observed in bulk GaAs at room temperature, in a recent pump-probe experi- 
ment with circularly polarized light. Hence our results may be useful in the microscopic modelling 
of density-dependent optical nonlinearities in spin-polarized semiconductors. 

PACS numbers: 71.10.-w, 72.25.Fe, 78.20.Ci 
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I. INTRODUCTION 

As the nascent field of spintronics^« is merged with op- 
toelectronics, it becomes increasingly important to un- 
derstand the physics of the optical properties of spin- 
polarized semiconductors, that are nowadays studied by 
circularly polarized pump-probe experiments.^ This de- 
gree of freedom of the spin, in the optical properties of 
photo-excited semiconductors, is mostly ignored in the 
literature. For spin-unpolarized photo-excited carriers in 
a pure semiconductor, it is known that the Coulomb at- 
traction between the electron and hole created by a pho- 
ton, causes excitonic resonances and enhances its absorp- 
tion;^ As the background electron-hole plasma density is 
increased, the long-ranged Coulomb interaction between 
this photo-excited electron-hole pair and the plasma can 
no longer be neglected. Semiconductor Bloch equations 
(SBE) fully describe this interacting plasma and its time- 
dependent dynamics. SBE may be solved in two dif- 
ferent time-scales: the coherent regime (where the in- 
duced dipole moment follows the optical field without 
de-phasing), or the quasi-equilibrium regime. In many 
applications, for example, semiconductor laser diodes, 
it may be assumed that the electrons and holes are 
thermalized within their respective bands."* This quasi- 
equilibrium approximation simplifies the original coupled 
SBE. The resulting microscopic theory for the effect of a 
spin-unpolarized interacting plasma (arbitrary density) 
on absorption, agrees well with experiment.^ Here the 
linear optical susceptibility is got by solving the Bethe- 
Salpeter equation in the (quasistatic) screened ladder 
approximation.^'^ In this formalism, the many-body ef- 
fects of the plasma on absorption are the screening of 
the Coulomb enhancement, phase-space filling (PSF) by 
the carriers and bandgap renormalization (BGR). These 
effects can be viewed as density-dependent optical non- 
linearities caused by the quasi-equilibrium plasmaj^ A 
two-band model is used for the conduction and valence 



bands by either ignoring the light hole [IK] band, or 
lumping it with the heavy hole {hh) band via an effective 
valence band density-of-statesi^ 

Photo-excited carriers may be spin-polarized by the 
optical orientation technique.^ It is based on the selec- 
tion rules for transitions induced by circularly polarized 
light, from both hh and Ih bands into the conduction 
band. A three-band SBE to model intervalence band 
coherence of quantum wells, under circularly polarized 
photo-excitation was formulatedji^ that included hh-lh 
band coupling. For modelling optical response in spin- 
polarized bulk and quantum well semiconductors, a very 
general six-band SBE has been recently framed and ap- 
plied (after neglecting some terms) to give a microscopic 
description of spin dynamics fii^i^ii^ As was done for 
the (spin-unpolarized) semiconductor laser^ii it maybe 
desirable to solve such model spin-SBE in the quasi- 
equilibrium regime, for a microscopic understanding of 
a new type of spin optoelectronic device, the spin verti- 
cal cavity surface emitting laser (spin VCSEL)ii^ 

In the quasi-equilibrium regime, using the pump-probe 
technique Nemec et al.l^ recently studied absorption 
spectra in spin-polarized bulk GaAs at room tempera- 
ture. They observed a spectral crossover in the difference 
in absorption between right ((T+) and left (ct^) circularly 
polarized light. This circular dichroism experienced by 
the probe, is due to the electronic spin-polarization ex- 
cited by an earlier pump pulse. For a microscopic de- 
scription of this experiment, we present the spin- modified 
solution of the Bethe-Salpeter equation, extended to in- 
clude the Ih band. Our approach is equivalent to solving 
the full spin-SBE in the quasi-equilibrium regime, but ne- 
glecting the terms corresponding to hh-lh coupling, spin- 
splitting of the single particle states and electron-hole ex- 
change interaction. The last two terms are important for 
spin relaxation processes, but not for optical transitions. 
hh-lh coupling was neglected because it is important only 
for inter-valence band processes. Numerical simulations 
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based on this framework, showed the spectral crossover 
and were in reasonable agreement with the experimental 
spin-dependent, density-dependent and spectral trends. 
The Bethe-Salpeter equation may therefore be useful in 
modelling spin-dependent many-body effects in semicon- 
ductors in the quasi-equilibrium regime. 



II. METHOD 

We first modify the solution of the Bethe-Salpeter 
equation, for spin-polarization (^) and inclusion of Ih. 
Then the method used to compare with the time- 
dependent experimental data^^ is described. 



Spin-polarized optical susceptibility (including 
light holes) 




FIG. 1: Selection rules in schematic band structure of bulk 
GaAs at the center of the Brillouin zone, with Wj > n|. In- 
cident light is assumed to be propagating in the z direction. 
Bands are labelled with their mj indices, where mj is the 
component of total angular momentum J along z direction. 
Solid and dashed arrows, with the relative transition rates in- 
dicated at their base, correspond to a'^ and a~ helicities of the 
probe. J, and j electron bands are labelled with mj = — 1/2 
and 1/2 respectively. The spin-unpolarized valence bands are 
labelled with mj = ±3/2 (hh) and ±1/2 (Ih). PSF at a non- 
zero temperature is shown by the gradation in thickness of the 
bands (the small PSF of the Ih band has not been shown). 
The greater BGR of the J, band causes its bandedge to be 
lower compared to the | band. 

At t — both electron and hole spins are created by 
the right circularly polarized pump pulse in the experi- 
ment; however hole spins relax within t ^ 100 fs leaving 
behind a spin polarization only from electrons. The se- 
lection rules favor creating spin-down (J,) electrons from 
the hh band three times as much as spin-up (t) electrons 



from the Ih band (J, and t electrons have their spin op- 
posite and along the propagation direction of the pump, 
respectively). Moreover the joint density-of-states effec- 
tive mass for hh transitions m^^ hh is nearly twice rrir^ ih 
for Ih transitions. These two factors cause = 6n|, as 
pointed out in Ref. [TgI (n^ — n, the plasma density). 
So, at t — the spin polarization ^max = 5/7, where ^ is 
defined as 



Later the same selection rules involving the hh and Ih 
bands, determine the absorption of the right ((7+) or left 
(cr") circularly polarized probe (Fig.[T]). Therefore, here 
the absorption depends not just on the total density of 
electrons and holes (as described in Refs. for 
^ = 0) but individually on n^, n|, rihh and nih- 

PSF reduces absorption since lesser number of states 
are made available for optical transitions. BGR in- 
creases absorption since the transitions take place at 
larger wavevector, as the bandgap is narrowed. PSF 
and BGR, which are density-dependent, are different for 
transitions into the | and "f electron bands (Fig. [1]) . As 
will be seen later, these transitions also have different 
Coulomb enhancement because the enhancement itself 
depends on PSF and BGR. All this affects the suscep- 
tibility x™{'-^) pertaining to optical transitions from the 
valence {v = hh, Ih) to the conduction (c =1, t) bands. 
The susceptibility x^(w) for can be written as a sum 
of two such transitions (Fig. [1]), 



X+(c.) = x^''"^(c.)±x^''''(c.) 
X-{u;) = x^''"^(c.)±x^'"'(c.). 



(la) 
(lb) 



is sensitive to the J, band, whereas x~ is weighted 
towards the t band, because the transition from the hh 
band is favored over that from the Ih band. Therefore, 
in general a circular dichroism results from a non-zero ^. 
The complex optical dielectric function e^(a;) is got from 
X*(cj) using 



±4^X^(w), 



(2) 



from which the absorption a^{uj) and circular dichroism 
Aa^Lu) = a'^^Lu) — a~{uj) were obtained. Also, since lin- 
early polarized light can be written in terms of and 
a~ , we obtain 



X 



X 



(3) 



where x*^ is the susceptibility of linearly polarized light, 
ft may be verified that if ^ = 0, x° = = X~ ■ 

The optical susceptibility x'^''(w) of a particular tran- 
sition was related to the microscopic susceptibility 
X™(fc,a;)by 



X™(^) 



(4) 
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The sum over wavevector does not include the spin de- 
generacy. The interband dipole matrix element for circu- 
larly polarized light, d™{k), reflected the selection rules, 



where Ago is the spin-orbit splitting of the valence band 
at the center of the Brillouin zone. 



d™(fc) 



ehPo 



moEg { 1 
'1 



ehPo 



3 moE„ ( 1 



[v = hh) (5a) 



[v = Ih), (5b) 



where mo is the free electron mass, to,-. „ = (I/toe -|- 
l/my)"^ and Eg is the bandgap. The momentum matrix 
element Pqj in an eight-band modcl^^ is given by 

Po^-^(^-l)g4|^, (6) 
2 \m,e J SEg + 2Aso 

I 



Interacting optical spectrum 



We repeat the steps for getting x'^'"{uj) outlined in 
Ref..l8„ but with modifications for Ih and spin. x™(fc, w), 
which causes the dichroism, is got by numerically solv- 
ing the effective Bethe-Salpeter equation describing the 
repetitive electron-hole scattering (ladder approxima- 
tion): 



[Uu; -Eg- A£;7 +*7-(fc,c.)] x""(fc,^) - -(1 - fc{k) - U{k)) 

I 



d^^(fc) + ^i4(|fc-fc'|)x™(fc'. 



(7) 



The screened Coulomb potential Vs{q) derived in the 
random phase approximation, and simplified in the qua- 
sistatic single plasmon-pole approximation is given by 



Vs{q) = 



1 47re2 



£3 gpg2 

47re^ 



1 - 



1 + 4 + f±^) 



47re^ dnj 
eo ^ dfij ' 



(8) 

(9) 
(10) 



where fij is the chemical potential, LOpi and k are the 
3d plasma frequency and wave number respectively and 
Ug simulates the electron-pair continuum. Note that Upi 
and K have been made ^-dependent. In Eq. (O, the factor 
{1 — fc{k) — fy{k)) causes PSF. The Fermi functions fj{k) 
that described the distribution of the electrons and holes 



fjik) 



exp[/?(i?,(fc)-^,)] + l 



(11) 



where Ej{k) = h^k^/{2mj) and (3 = l/iksT). We ig- 
nored the small spin-dependence of the electron massif 
(i.e. TOj — mij = nie)- The BGR, Ai?™, for an optical 
transition between the valence band and the conduction 
band is 



AE^^k) = edk) + e^k). 



(12) 



The self-energy ej{k) of the j quasiparticlc, in the qua- 
sistatic approximation is 

e.-W :^-Y.Vsi\k-k'\)f,{k') + l[Vsir = 0)-V{r = 0)], 

k 



(13) 

where the first term is the 'screened-exchange' and the 
second term is the 'Coulomb- hole'. The Coulomb-hole 
term is the same for each of the j quasiparticles. In 



Eq. 



which simulates the electron-pair continuum. 



as well as the temperature-dependent damping 7'^"" (k,ijj) 
and bandgap Eg in Eq. ([7]) were got from Refs. [l^ and[^. 
The chemical potentials fij, assuming that the electrons 
and holes were in quasi-equilibrium, was got using a form 
of the Aguilera-Navarro approximation.— The hh and Ih 
have the same chemical potential because they are in 
equilibrium with each other. 
Rearranging Eq. ([7]) we got 



X™(fc,t^)=X§"(fc,c.) 



where 

X'o^'ik^uj) 



(14) 



d--ik){l-Uik)-U{k)) 
hw- Eg- AE^" + i Y^{k, uj) ' 



(15) 



To get the correct crossover between gain and absorp- 
tion with 7™(A:, lu), Xo^{k, w) was described by a spectral 
representatiouji^ By defining 



X™(fc,c.)=r-(A:,c.)xS'^(fc,c.) 



(16) 
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and substituting it into Eq. (|T4)) . we obtained for the 
vertex function T™(k,u!), 



tance D 



r™(fc,w) = 1 + 



k' 

(17) 

Vs{\k — fe'l) has been replaced in Eq. ([TT]) by its angle- 
averaged value Vs{k,k'), because we assumed that only 
s-wave scattering contributed to the optical transitions. 
We used a matrix of approximately 200 x 200 Gauss- 
Legendre quadrature points to represent the vertex inte- 
gral equation (flT)) . The diagonal singularity of the matrix 
was regularized by the compensation technique/ before 
it was inverted to give the solution for r'^"(fc, uj). The op- 
tical susceptibility x'^"(a;) of a particular transition was 
obtained from 



^^d-(fc)r™(fc,c.)xS^(A:,c.), 



L3 



(18) 



where we substituted Eq. (HH) in Eq. g]). 



2. Noninteracting optical spectrum 

If r^"(fc,a;) is neglected in Eq. we obtain the 

'noninteracting' susceptibility, which differs from the 
truly noninteracting susceptibility due to the BGR term 
present in x[5''(fc,w) [c.f. Eq. 1^]. r™(fc,w), which de- 
pends on PSF and BGR [c.f. Eq. ^ and ([15])], causes 
the excitonic resonances and the Coulomb enhancement. 
This is because r™{k,uj) expresses the influence of mul- 
tiple electron-hole scattering (resulting from their attrac- 
tive interaction), on the susceptibility. 

It is difficult to separate PSF and BGR in the in- 
teracting Aa(w) due to the presence of the r™{k,uj) 
term. Therefore we used the noninteracting x™{'-^) to 
study how competition between PSF and BGR influences 
Aa(w). 

In our calculated spectra, we took the spectral rep- 
resentation for only the imaginary part of xg^(A:, w) 
Also, the fc-dependent BGR was taken as a rigid shift at 

1/3 

kp = (Stt^ti) . These simplifications are not expected 
to significantly affect the results. The material parame- 
ters used were: exciton Rydberg Eo — 4.2 meV, exciton 
Bohr radius ao = 125 A, electron mass me — 0.0665 mo, 
hh mass rrihh = 0.457 mp, Ih mass mih = 0.08 mg, 
Aso = 0.341 eV eo 13.71, e^o = 10.9 and T ^ 295 K. 



B. Comparison with experiment 

The probe spectral width was accounted for by adding 
its 'half-width at half-maximum' value of 15 meV to 
7™(fc,Ci;) in the calculation. The experimental data is 
mostly in terms of the normalized differential transmit- 



D = 



\ T I 



(19) 



Here {AT/T)^ = {T^-T)/T, where T± is the transmis- 
sion of probe cr^ after the sample is pumped with cr+ . T 
is the unpumped transmission through the bulk sample 
of thickness w ~ 1/im. If the change from the unpumped 
absorption Aa^ is such that Aa^ <^ 1/w ^ lO'^cm^^, 
we can write D as 



D 



Aa' 



Aa' 



2ao 



2ao 



(20) 
(21) 



The unpumped absorption aQ{uj) could be calculated 
from x^{^) with ^ = and the background doping den- 
sity of the sample 



15 



no 



10^ cm""*. The sign of D is 
opposite to that of Aa because screening by the pumped 
carriers usually reduces the absorption i.e. a'^ + a~ — 2ao 
is a negative quantity. Furthermore, D is not affected by 
an overall density-independent scaling factor C{uj) that 
could multiply the calculated absorption, since Z? is a 
normalized quantity. Usually C{uj) is needed to match 
the calculated absorption to the experimental value of 
absorption in a 'pure', unexcited sample.^ 

The ^-dependent calculations were compared with the 
time-dependent experimental data by using the following 
equation describing spin relaxation. 



C(^) — ^max 6 



~2t/T 



(22) 



where 



5/7 and the spin relaxation time*^ r. 



130 ps. We assumed a temporally constant plasma den- 
sity n since f ^ r,. « 1 ns, where is the carrier recom- 
bination time. 



III. RESULTS AND DISCUSSION 

We show that our calculation captures all the trends 
in the experimental data, in probed energy hv, spin po- 
larization ^ and pumped density n. To explain how PSF 
and BGR cause a (^-independent) crossover in Aa, we 
first discuss results for noninteracting Aa. After show- 
ing the effect of the Coulomb interaction on Aa, we fi- 
nally compare the results for interacting Aa with the 
data. The agreement with the data comes directly from 
our calculation, without requiring us to adjust ^, n or 
sample-dependent broadening. 

We did two checks on the spin-polarized calculation. 
Putting ^ = gave back the unpolarized absorption 
(Fig-Uta)) as calculated by Ref. llSl for different plasma 
densities. At lower densities, with ^ = 0, we also got 
the broadened Elliott's formula^ by removing the k- 
dependence of '-f™{k,uj) and d™{k). The discrepancy 
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with Elliott's formula was « 10"'^ for photon energies 
that exceeded the bandgap. 

Inclusion of Ih serves only to enhance the unpolar- 
ized absorption by 1/6, 1/3 from the matrix element and 
1/2 from the density-of-states (compare dashed and solid 
curves for n — noin Fig.[2ja)). Treating Ih independently 
(without hh-lh coupling) did not lead to noticeable arti- 
facts in the calculated absorption. 



is caused by competition between PSF and BGR. Within 
a transition (either hh or Ih), PSF dominated at lower 
energies whereas BGR (of the opposite sign) dominated 
at higher energies. Moreover, PSF and BGR from Ih 
transitions were opposite in sign (due to the selection 
rules in Eq. ([1])) compared to those from hh transitions 
(Fig. m^b)). The Aa crossover shifted to higher energies 
as n was increased, because PSF became more impor- 
tanO^ (Fig.IHc)). 



A. Noninteracting optical spectra 

The microscopic noninteracting calculation agrees with 
the earlier simplified explanatioi*^^ that the Aa crossover 
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FIG. 2: (a) Spin-unpolarized interacting a, with background 
no = 10^^ cm^'^ and photo-excited plasma densities Fo — 
1.3 X 10^^ cm^'^ and SFq. The dashed curve differs from the 
solid curve with n — no, by the neglect of Ih. (b) Nonin- 
teracting Aa decomposed into hh and Ih transitions, each of 
which has contributions from PSF and BGR. Within a transi- 
tion, PSF and BGR oppose each other. Between transitions, 
PSF and BGR from hh oppose those from Ih. (c) Density- 
dependence of total noninteracting Aa. The curve for n = Fq 
can be got by summing together the curves in (b). 
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( hv - E J / E„ 



FIG. 3: (a) ^-dependence of noninteracting Aa. (b) PSF 
and BGR components of noninteracting a'^ at the crossover 
energy vs. (neglecting Ih band). The ordinate is normal- 
ized to the value of a^ at n| — n/2, whereas the abscissa 
is normalized to n. The opposite n|-dependence of PSF and 
BGR add up to give total a^ (scaled by 0.5) that is almost 
n^-independent. The ^-independence of Aa at the crossover 
energy in (a) may be got from the difference of the value of to- 
tal a^ at the values shown by vertical solid lines (^ — 5/7), 
dashed lines (^ = 0.5) and dot-dashed lines (^ = 0.16). 



In contrast, the Aa crossover energy was independent 
of C at a fixed n (Fig. due to the PSF and BGR 

interplay. A similar ^-independent crossover occurred at 
higher energies, if hh transitions alone were considered. 
The weaker Ih transitions did not result in a crossover 
in Aa in the energy interval shown, as can be verified 
by summing the dashed curves in Fig. [IJb). Since ^ de- 
termines ni relative to n|, we study the cause of the 
^-independent crossover by plotting a"*" (inducing transi- 
tions only between the hh and | bands after neglecting 
the Ih transitions) vs. (Fig. [3]Jb)). Analogous results 
were got by plotting a~ vs. n|. The n|-independence of 
total a"*" at the crossover energy, indicates that opposing 
trends between PSF and BGR for each helicity of light, 
caused the Aa crossover to be independent of i^. Only 
the electronic part of the Ih transitions contributed to 
Aa (holes are unpolarized), hence the ^-independent na- 
ture of the crossover was unaffected in the full Aa which 
included Ih transitions (Fig.[3Ja)). 
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B. Interacting optical spectra 
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FIG. 4: (a) Density-dependence of interacting Aq. Tlie cor- 
responding noninteracting Aa sliown, is identical to Fig.[2jc). 
(b) ^-dependence of interacting Aa. 

The Coulomb interaction increased the magnitude of 
Aa and shifted its crossover energy, but did not affect the 
trends in hv, n and ^. The interaction caused a sharp 
peak near the bandgap (due to excitonic effects) in the in- 
teracting Aa, for n — Fq, and also enhanced it at higher 
energies compared to the noninteracting Aa (compare 
solid and dashed curves for n = Fo in Fig. [D^a)). At 
n = 8Fq, the peak near the bandgap in the interacting 
Aa was not as pronounced (due to the almost complete 
excitonic ionization) compared to the noninteracting case 
(solid and dashed curves for n = SFq in Fig.lHa)). But at 
higher energies, they show a persisting Coulomb enhance- 
ment even sd, n — 8Fq. This is because the screening by 
the plasma is not very effective at high energies^ (solid 
curves for n — uq and 8Fq in Fig. [U^a) approach each 
other at high energies). However the ^-independence of 
the Aa crossover energy was preserved, despite Coulomb 
interactions (Fig. Sfb)). The cause of the shift of the 
crossover to lower energies when we compare the inter- 
acting with the noninteracting curves is unclear. 



C. Comparison with experiment 

We compare our calculation with the experiment by ac- 
counting for probe width in the broadening (Fig. [Sla)). 
We think that spectral averaging due to this additional 
broadening causes the crossover to shift from {hv — 
Eg)/Eo = 10.5 (solid curve for n = Fq in Fig. IJa)) 
to {hv — Eg)/Eo = 15.3. Because D is plotted instead of 
Aa in Fig. ^a,) , there is an overall flip of sign compared 
to Fig. H 



Our model reproduces the experimental observations 
that D changes sign either as A is varied at a fixed n 
(Fig.[5fa)) or n is varied at a fixed A (Fig.[5]^b)). Further- 
more, the expected ^-independence of the crossover en- 
ergy (Fig.|4l^b)) is indeed shown by the data (flat dashed 
curve in Fig. [S]). In Figs.jHJb) and [S] the calculation used 
A ~ 816 nm and A — 834 nni instead of the actual probe 
A = 775 nm and A — 800 nm respectively. The tendency 
of D to change sign at a fixed probe wavelength or en- 
ergy (Fig.EJb)) as n is increased, occurred over a certain 
range of energies. This tendency is also shown by Aa 
over 10.5 < {hv — Eg) / Eq < 23.1 (solid curves for n — Fq 
and 8Fo in Fig. [4{a)). Outside this range \D\ increased, 
without flipping its sign, with n. We cannot obtain the 
carrier thermalization part of the experimental spectra 
{t < Ips) in Fig. [S] because it has not been accounted for 
in the calculation. 

The maximum value of D due to PSF is 0.35 (instead 
of 0.25 as stated in Ref. T^. This is because, by assum- 
ing Aa ~ —n in Eq. pp)) . we obtain D — and as 
noted in Section [TTl ^,nax = 5/7. The numerically calcu- 
lated D supported the above reasoning: in Figs. and [51 
D ^ 0.35 at lower energies. However near the bandgap, 
perhaps the excitonic nonlinearity slightly increased the 
value of D (Fig.lSJa)). \D\ significantly exceeds 0.35 only 
at higher energies where PSF is expected to be less im- 
portant (since D is negative) . The quantitative mismatch 
with experiment, at higher energies may be due to the 
neglect of dynamic screening and band non-parabolicity. 
We adjusted within a factor of 2 the values of v^ and the 
numerical constant a (describing ^{k,uj))^ but found 
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FIG. 5: Comparison with experiment, (a) Spectral D. The 
data points are got from Figs. 4(b) and 6(b) of Ref. [Tsl , at 
t = 7 ps. (b) Density-dependence of D at a fixed probe energy. 
The experimental points are for n = Fq, 2Fo, 4Fo and 8Fo, 
again at t = 7 ps (from Fig 7 of Ref.lTsh. The calculation used 
A — 816 nm whereas experimentally A — 775 nm. Vertical 
bars indicate the experimental noise. 
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that the crossover energy changed by ^ 1%. Experimen- 
tal uncertainties in the photo-excited carrier density n 
and the induced spin-polarization ^ as well as sample- 
dependent broadening, none of which were adjusted for, 
are also expected to affect the mismatch. 
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IV. CONCLUSION 

We have provided a method to calculate spin-polarized 
many-body effects in the room temperature absorption 
spectra of bulk GaAs. This was done by modifying the 
existing microscopic theory for absorption, that is known 
to match experiment over a wide range of plasma den- 
sities. Light hole contributions were also included. The 
agreement with results of a recent circularly polarized 
pump-probe experiment came directly from our calcu- 
lation, without requiring us to adjust spin-polarization, 
plasma density or sample-dependent broadening. We 
find that the Aa crossover and the experimental spin- 
dependent, density-dependent and spectral features are 
reproduced, thus validating the use of our model to un- 
derstand circularly polarized pump-probe experiments in 
III-V semiconductors, in the quasi-equilibrium regime. 
This also opens up the possibility that the Bethe-Salpeter 
equation may be used to theoretically describe spin- 
dependent many-body nonlinearities in the operation of 
the spin VCSEL and perhaps make predictions about its 
performance. 
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FIG. 6: Comparison with experiment (dot-dashed lines) re- 
garding evolution of D at the probe A = 775 nm, 800 nm (834 
nm in the calculation) and 865 nm. The experimental data 
are got from Fig. 8 of Ref . [isi. Pumped densities are 1.3 x lO^'^ 
cm"^, 6 X lO^*' cm~^ and 2 x 10^*^ cm"^ respectively.— Vertical 
bars indicate the experimental noise. 
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